Topological and geometric decomposition of nematic textures 

Simon Copar 

Faculty of Mathematics and Physics, University of Ljubljana, 
Jadranska 19, 1000 Ljubljana, Slovenia 

Slobodan Zumer 
Faculty of Mathematics and Physics, University of Ljubljana, 

Jadranska 19, 1000 Ljubljana, Slovenia and 
Jozef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia 
(Dated: January 23, 2012) 

Abstract 

Directional media, such as nematic liquid crystals and ferromagnets, are characterized by their 
topologically stabilized defects in directional order. In nematics, boundary conditions and surface- 
treated inclusions often create complex structures, which are difficult to classify. Topological charge 
of point defects in nematics has ambiguously defined sign and its additivity cannot be ensured when 
defects are observed separately. We demonstrate how the topological charge of complex defect 
structures can be determined by identifying and counting parts of the texture that satisfy simple 
geometric rules. We introduce a parameter called the defect rank and show that it corresponds to 
what is intuitively perceived as a point charge based on the properties of the director field. Finally, 
we discuss the role of free energy constraints in validity of the classification with the defect rank. 

PACS numbers: 02.40.-k,61.30.Jf,61.30.Dk 
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I. INTRODUCTION 



Various states of matter we know today are characterized by interactions and symmetry 
of their microscopic building blocks. Local order is reflected in mesoscopic quantities, such 
as the phase of electromagnetic waves, crystalline order, nematic director, magnetization, 
dielectric tensor and others. We call these quantities order parameters. A prominent example 
of partially ordered matter are liquid crystals, which are characterized by orientational 
order and liquid-like positional freedom of the molecules. Nematic liquid crystals exhibit 
topologically prescribed defects, which in combination with surface-treated colloidal particles 
nanoparticles or general confinement 0, enable creation of complex materials with 



tunable optical response [l]. 

The nematic order can be represented with an order parameter tensor that contains both 
the magnitude of the ordering and the directional information In studies of defect topol- 
ogy, the principal eigenvector of this tensor, the director, is a common choice to represent 
the average direction of molecular alignment. Vectors with unit size lie on a sphere {S 2 ), 
which is a topologically nontrivial space. Consequently, boundary conditions can exist that 
cannot be smoothly interpolated in bulk without introducing defects 6|, |7| . In contrast to 
Heisenberg ferromagnets, which have the ground state manifold with exact topology of a 
sphere, nematic liquid crystals have an additional rule that opposite vectors represent the 
same state, as the molecules have no back-front distinction [8|. One consequence of this 
property is that beside the point defects, nematics also allow line defects. Beside that, the 
sign of the vector field representing the director is ambiguous and therefore the same holds 
for the invariants of the field which are odd in the director, such as the topological charge 



of point defects 

The growing interest for controlled formation of complex defect structures in constrained 
nematic liquid crystals [l-3] is increasing the demand for visualization and simple classifica- 
tion of complex nematic director fields. This stimulated us to derive a natural expression for 
the topological charge that does not involve integration and reflects its discrete nature. We 
accomplish this by decomposing the director field into a hierarchy of geometric primitives, 
which enables intuitive understanding of the topological charge based on the qualitative 
features of the director field. We extend the classification from strictly topological to a more 
general one that has the advantage of differentiating topologically equivalent states with 



different energies. We introduce a quantity called the defect rank and draw parallels with 
the topological charge and naming conventions for point defects. 



II. PATCHES, BOUNDARIES AND GRAINS 

If the vector heads are artificially prescribed, the director field n becomes a proper vector 
field. Volumes of interest, such as those that include defects, can then be classified using the 
topological charge, which enumerates the defects with the elements of the second homotopic 
group p, |8j. As topological properties are invariant to the specific shape of the surface that 
encloses the test volume, a spherical volume is assumed from here on without the loss of 
generality. For a sphere, parametrized by the spherical angles 9 and 0, the Gauss integral 
for the topological charge [3, B-tlll reduces to the following expression, 
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This expression involves integration over all the details of the director field n, which can be 
eliminated, if the director n on the enclosing sphere (the surface texture) is decomposed into 
geometric primitives, over which the integration can be performed in advance. According 
to the component of n along the surface normal, the sphere can be split into patches of 
outgoing vector field and patches of ingoing vector field. These patches are separated by a 
disjoint set of closed loops (boundaries) where the director lies tangentially to the enclosing 
sphere. On the boundaries, the director still has one degree of freedom - its direction in 
the tangent plane. With the same approach as for the patches, the boundary can be split 
into parts with a constant orientation along the normal to the boundary curve. These are 
separated by points (grains), where n is directed tangentially along the boundary (Fig. [Ik). 
Degenerate cases, such as textures with touching boundaries or wide areas where director 
is parallel to the surface, can always be regularized by an infinitesimal perturbation of the 
director field. 

The parametrization of the sphere with the polar angle 9 and the azimuth generates a 
local coordinate frame (f, 9, 0), where 9 and lie parallel to the surface (Fig. UK). Within this 
coordinate frame, we introduce local spherical coordinates with angles a and {3 to express 
the director n, 

n = t(f cos a + 9 sin a cos + sin a sin (3) (2) 




FIG. 1. Geometric decomposition of a surface texture on a sphere, (a) Any texture can be smoothly 
deformed to make the director perpendicular to the sphere everywhere except in narrow transition 
zones. This decomposes the surface into patches with perpendicular director, boundary loops with 
parallel director and point-like grains with the director rotating in-plane for half a turn. The inset 
shows the continuous texture of a grain, (b) Grains contribute g = ±1/2 to the topological charge, 
depending on the sense of the director rotation. The table shows four possibilities, together with 
corresponding (t, s, w) parameters that define the director's orientation in patches, boundaries and 
grains, respectively, (c) Merging and dividing of patches by continuous transformations preserves 
the topology of the texture. The change in the number of patches is compensated by spontaneous 
creation of grains. 

where t = ±1 simplifies the calculation for director fields with reversed vectors. Without the 
loss of generality, each angle can vary in only one of the directions: a = a(8) and (5 = /?(</>). 

To calculate the contributions of patches, boundaries and grains to the topological charge 
(Eq. [T]), a model for a(8) and /3(<f>) is needed. On every patch, the director n can be combed 
perpendicularly to the surface without changing the topology, because the condition for 
the normal component to have a constant sign locally restricts the ground states to a half- 
sphere, which is topologically trivial (contractible to a point). This is achieved by setting 
a = (3 = so that the director points directly outwards if t = 1 and inwards if t = — 1. For 
a single patch, shaped like a spherical cap bounded by a circle of constant latitude at 9 = 9 
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(Fig. the contribution to the topological charge follows, 

q ' = T«l I H^ x ¥) d( ' d * = 5 ((1 - C08( '» ) ' (3) 

This contribution is proportional to the fraction of the sphere covered by the patch and its 
sign depends on whether the director points inwards or outwards from the surface, which is 
given by the signature t = ±1. 

On the boundary, the director makes a smooth turn from outward to inward direction. 
A boundary without grains is modeled by setting (3 = and choosing a to do a half-turn 
within a narrow transition zone between the boundary and the patch, 

°<«=-K i+ ^)- (4) 

7] is the half-width of the transition zone around the boundary, and the boundary signature 
s = ±1 selects the in-plane orientation of the director at the boundary, n = st9 (Fig. [Th). 
To simplify the calculation, the width rj of the transition zone is limited to zero, 
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The result does not depend on s, therefore the orientation of the director on the boundary 
does not influence the topological charge. 

The contribution of a single grain is obtained in a similar fashion. The reverse of the 
in-plane orientation of the director is smoothly interpolated by (3(<f)) = iw7r/2(l + (0 — <j>o)/£), 
with £ being the half width of the transition in direction (Fig[Th,inset). Parameter w = ±1 
governs in which way the director points at the grain, n = (tsw)<fr. A limit £ — > yields the 
contribution of the grain. 
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Consider a texture, divided by N boundaries into N + 1 patches of arbitrary shapes. 
Every boundary can have an even number of grains, since the continuity requires an even 
number of director flips on a closed circuit. The topological charge q of the entire texture 
then splits into a sum of elementary contributions (Eqs. 

N N M 

q = ^2 ( ipi + ^2 ( ibi + Yl 9i- ( 7 ) 




boundaries 



grains 



FIG. 2. Boundaries divide the texture into patches with arbitrary positioning, shape and hierarchi- 
cal nesting order. If the texture is instead viewed as a superposition of layers, the boundaries can 
be counted together with patches. In such decomposition, each patch contributes ±1 to the topo- 
logical charge, depending on whether its direction is oriented inwards or outwards, here represented 
with alternating colors. 



Instead of adding patch contributions q P i, which are proportional to the areas of the patches, 
multiplied by the signature t (Eq. [3]), we start with a uniform background sphere and pro- 
gressively build the desired configuration. For convenience, we choose the background with 
the director pointing outwards, t — 1. To recover the desired configuration, N patches must 
be added. Each patch contribution is counted twice, once to cancel out the background and 
once to account for its own contribution. This way, all the patches are shaped as spherical 
caps (Fig. [2]) and their contributions can be grouped with their respective boundaries, 

N N 

q + ^{2q pi + q u ) = 1 + ^ U- (8) 



N N 
i=l 



i=0 i=l i=l i=l 

The background patch q has the topological charge of a radial hedgehog, q 



X/2 



where x is the Euler characteristic of the surface |12l.ll3|. The dependence on the polar angle 
9 has disappeared, since the shape and position of the boundaries are not topologically 
relevant. Any number of boundaries can be nested one inside another, which results in 
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patches with alternating t being added to the sphere (Fig. [2]). 

The grain contribution g = tw/2 can be associated with geometric features of the texture. 
The director n rotates around the surface normal as we move along the boundary. There 
are two possible directions in which we can trace the boundary, but they are invariantly 
distinguishable, as there are patches of opposite signatures on each side. For a choice of 
the direction that has the positive signature patch on its left side, the grain contribution is 
+ 1/2 if the director rotates anticlockwise and —1/2 otherwise (Fig. [lb). 

The charge decomposition in the form 

N M 

<Z =! + !> + !> (9) 

i i 

applies not only to the nematic director, but to any medium that can be described with a unit 
vector order parameter. It works for spherical enclosing surfaces, with possible generalization 
to surfaces of arbitrary genus. We assumed the vector field on the enclosing surface is non- 
singular. Singular objects, like boojums and other surface defects, contribute additional 



terms to the expression 
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Patches, boundaries and grains are two-, one- and zero-dimensional geometric features of 
the surface texture. This decomposition is similar to the notion of domain walls and defect 
lines in ferromagnetic systems jla \\\ , except in our case we observe 2D texture on a closed 
surface instead of bulk 3D field. Our boundaries are similar to cross sections of magnetic 
Neel walls and the grains correspond to cross sections of Bloch lines. However, in magnetic 
systems, the domain walls and lines naturally prefer narrow transition zones 16[, whereas 
in nematics narrow transition zones are merely a computational tool. 



III. FREE ENERGY CONSTRAINTS 



One of the reasons for calculating the topological invariants is to detect similarity between 
director fields. Take two surfaces that enclose closed volumes in the bulk. If we cut out the 
volumes and exchange them, we can only reconnect them smoothly with the ambient director 
field if they are topologically equivalent. Topological charge tests for compatibility of bulk 
volumes of the medium judging only by the director on their surfaces, so the charge can be 
assigned to the surface texture itself. Compatibility of two surface textures neither implies 
additivity of the topological charge hidden inside nor the invariance to the choice of the 

7 



enclosing surface. We are free to add further restrictions and thus refine the classification. 

Despite the topological equivalence of two volumes, the transition between them may not 
be possible due to the energy constraints, which depend on the medium and the choice of 
the enclosing surfaces. The structures can be deemed equivalent if one can be transformed 
smoothly into the other without deviating significantly from the local free energy minimum 
(if there is a transformation that avoids major free energy barriers). The deformation part of 
the free energy is related to the spatial derivatives of the director. A general transformation 
can both deform the shape of the medium and independently rotate the director locally. 
However, if the director is treated like it is pinned to the underlying coordinate space and 
rotates with it, the neighborhood of each point stays approximately the same, except for 
the anisotropic distortions, which can be considered small. By allowing only the subset of 
all smooth deformations that obey this rule, topologically equivalent states with different 
energies can be differentiated. In topology, such mappings are called push-forward mappings. 
This type of restriction is more natural for solids, where the underlying space corresponds 
to the crystal lattice. We can argue that two structures are similar if we could deform one 
into the other if they were elastic solids. 

The surface texture decomposition involves expressing the director in a coordinate frame 
aligned to the surface (Eq. |2J). Under a mapping that satisfies the restriction to push- forward 
transformations, the local coordinate system retains its relative position to the director and 
the decomposition does not change. The texture decomposition is therefore a property that 
classifies the surface textures into classes, where the representatives of the same class can be 
stitched together without significant change in free energy. The topological charge, which 
can be computed from the decomposition (Eq. [9]), is only part of the information included 
in the full decomposition. Other invariants, like the relative position of the boundaries and 
their boundary signatures s — ±1, which have no effect on the topology, are also good 
for classification of the qualitative shape of the texture. In the following section we show 
that these quasi-invariant s are relevant for nematics, where traditional topological charge is 
problematic due to the additional symmetries of the director. 
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IV. DEFECT RANK OF NEMATIC TEXTURES 



Topological charge (Eq. [T]) has an odd parity under inversion of the director n h-> — n. The 
additivity of the topological charge, which holds for true unit vector fields, may be violated 
because nematics allow line defects with half-integer winding number. These defects cannot 
be represented as a vector field without explicitly introducing a branch discontinuity where 
the vectors flip sign. In the presence of line defects, the topological charge of each individual 
point defect can be flipped just by bringing it across a branch discontinuity, without changing 
its physical shape p, |7|. The topological charge of a cluster of point defects cannot be 
determined from the individually measured charges of constituent point defects. Instead, 
only even or odd parity of the topological charge is properly defined and conserved [isl . 

Despite ambiguity, the signs of topological charge are usually assigned by convention. The 
charge +1 is assigned to all radial hedgehogs and —1 to hyperbolic hedgehogs, regardless the 
orientation of the arrows on vectors, even though there is a 3D rotation in director space that 
can transform one into the other In case the hedgehogs interact in the same medium, 

the continuity of the director field prescribes their relative signs of topological charge that 
may clash with this convention. The transformation from a hyperbolic to radial hedgehog 
involves a rotation of the director without simultaneous rotation of the underlying space. 
Such transformation does not satisfy our condition to push-forward transformations and the 
incompatibility of the surface textures means they cannot exist in the same surrounding 
director field without costly elastic deformations. The transition between the hedgehogs is 
thus restricted by an energy barrier, except in cases where the experiment is specifically 



designed to support the transition 



19] . 



We seek a geometric invariant that matches the conventional understanding of the hedge- 
hog charge and preserves even or odd parity, without introducing any inconsistencies not 
already present in the topological definition. In the following section, we define the defect 
rank of a point defect, q*, and compare it to the topological charge q. 

The definition of the defect rank q* can be made consistent with the conventions if we 
find a matching integer expression involving the surface texture properties. Earlier (Eq. H|) 
we introduced the boundary signature s = ±1 that describes how the director behaves on 
the boundary. The boundary signature s = —1 corresponds to the director stream lines that 
touch the boundary from the outside of the sphere (Fig. [3^) and s — 1 to the streamlines that 
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FIG. 3. The geometric meaning of the boundary signature s. (a,b) The boundaries are the 
curves where the streamlines of the director field touch the enclosing surface. The boundary 
signature s distinguishes whether the critical streamline (thicker line, red online) lies inside or 
outside the enclosing surface. It is not a topological invariant, but is invariant to the director sign 
reversal. The streamlines show the director in a radial cross section of the sphere, (c) Schematic 
representation of two surfaces with nematic textures in the cross section. Surface textures are 
compatible if their boundary signatures s match. The boundaries of opposite signatures cannot be 
interpolated smoothly, unless the director is rotated through the third dimension while the surface 
stays undeformed, which is forbidden in our classification rules. 

touch the boundary from the inside (Fig. [3b). This behavior is invariant to the inversion of 
vectors, therefore s is a good parameter for classification of nematic textures. Whether the 
streamlines that go through the boundary curve lie inside or outside the enclosing sphere 
can be tested by observing how the normal component of n along the direction of n. A 
local expression s = — sgn((n • V)(n • u)) can be constructed, which assigns the boundary 
signature to each point of the boundary and is also invariant to the inversion of n. 

Two surface textures are compatible if all the boundaries and grains are similar and 
arranged in a similar hierarchical pattern. To extract an invariant that resembles the point 
charge, we disregard most of the information and only focus on the boundary signatures 
s. Grains are the points on the boundary where s changes sign. The boundary signature 
s can be assigned to a boundary as a whole only if there are no grains present on the 
texture, otherwise it is a local property that is only constant between two grains. To find 
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a suitable expression for the defect rank, the classification must be restricted to textures 
without grains. Two such textures match if their boundaries have the same signature s, 
otherwise the stitching involves rotation of the director relative to the local coordinate 
frame (Fig. [3b). The topological charge (Eq. alternates between even and odd with the 
number of boundaries. This property is replicated by the defect rank, if we define it as 

N 

q * = i + Y^ Sl . (io) 

i 

If the surface has grains in its texture, the defect rank q* is not defined. The textures 
with grains can be understood as intermediate states through which the texture must pass 
in order to change its defect rank. For the defect rank to be a useful parameter, such 
intermediate states must have a higher energy, although the energy constraints are not 
absolutely prohibiting, like topological are. If the energy barrier is removed, which can be 
accomplished by external fields, boundary frustration or heating, the transition occurs and 
the texture changes Q, ljj. The entire texture on the enclosing surface can still be used 
for more precise comparison and testing for compatibility. The texture decomposition can 
be seen as a generalization of the traditional topological charge, as the latter is included in 
the former, as given by (Eq. |9]). 

The boundary signature s is not a topological invariant, but nevertheless has interesting 
properties which resemble conservation laws needed for effective classification. Two neigh- 
boring boundaries with opposite signatures can annihilate, as they represent a "fold" in the 
director field (Fig. [4b,). The definition of q* matches the convention for the radial hedgehog, 
which has no boundaries, (q* = 1), and hyperbolic hedgehog that has two boundaries with 
negative signature (q* = 1 + 2 ■ (— 1), Fig. lib). The singular case when the boundary is 



shrunk to a point corresponds to a point defect sitting on the enclosing surface [15]. Pushing 
additional point defect with rank s inside the enclosing surface adds a boundary of the same 
signature s to the surface, which suggests that the defect rank is conserved. 

Consider the case of two radial hedgehogs arranged into a vector field resembling that of 
an electric dipole. The hedgehogs must have the opposite orientations of arrows, as one acts 
as a source and the other as a sink. This structure has q = by strict topological definition 
(Fig. 0b). However, we observe that the texture on an enclosing sphere has one positive- 
signature boundary, which translates to q* = 2, the same as we get with naive counting 
of hedgehogs (Fig. 0H). Note that the "dipole" in this case is not the structure with one 
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FIG. 4. Point defect charge in nematics visualized with director field streamlines in cross section, 
(a) The signatures of boundary pairs, caused by director field bending are opposite and conserve 
the defect rank, (b) Hyperbolic hedgehog is identified by the defect rank of q* = — 1. (c) For the 
dipole configuration, strict topological definition assigns the opposite charges to the two hedgehogs, 
corresponding to the direction of the arrows and assigns the charge q = to the whole structure, 
(d) The defect rank treats the hedgehogs equally and amounts to q* = 2, the same as expected 
with naive counting of the hedgehogs. The arrows are not needed for determination of the defect 
rank. 

hyperbolic and one radial hedgehog, which is usually denoted by this name in nematics. 
The defect rank detects this and differentiates between them, while the topological charge 
assigns q = to both. 

The defect rank q* is sign-invariant and matches the existing conventions [k]] for assigning 
the charge to the structures in nematics. Despite similarities, this quantity is not the 
topological charge and only stays invariant under restricted set of transformations that 
preserve the texture on a chosen surface. The choice of the surface is no longer arbitrary: 
different surfaces may result in different q* even if they enclose the same topological charge. 
The difference appears when the contents of one surface do not fit into the other surface 
without deviating from the energy minimum. 

The simplicity of texture decomposition and defect rank can be demonstrated on the 
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FIG. 5. Two examples of nematic textures with marked patches (red and blue), boundaries (yellow) 
and grains (green), (a) Analytic model of a nematic point defect with q = —4. The texture has 5 
patches with t = — 1, which amounts to q = —4. The boundaries have signatures s = — 1, which 
gives the defect rank of q* = —4. (b) Model of a nematic point defect with q = 4. There are 3 
patches with t = —1 and 12 grains with g = + Together, this reproduces q = 4. However, the 
defect rank cannot be assigned due to the presence of grains. 



model director fields for point defects with known topological charge 10|, |20|, |21] . Visual- 
ization of the normal component of the director, n • u, shows that negative q solutions are 
realized as \q\ + 1 negative-signature patches on the enclosing sphere (Fig. [5^). For these 
solutions, the defect rank can also be determined and is equal to the topological charge. On 
the other hand, positive q solutions have q — 1 negative-signature patches and each boundary 
has four g = +~ grains, which gives the correct topological charge (Fig. |5b), but the defect 
rank cannot be determined due to the presence of grains. In nematics, the entire texture can 
be arbitrarily multiplied by —1, so the q and — q textures are topologically equal. However, 
the surface textures and their defect ranks are different, which implies that the two states 
are not interchangable without deviating from the energy minimum. Textures with high 
defect rank or topological charge are energetically expensive and as such hard to find in real 
systems. In an idealized calculation, the only stable point defects are those with q = ±1 



and in a cluster of many point defects, they prefer to form pairs with zero total charge 22 ]. 
Complex textures are most likely to be found locally as a part of a larger complex of line 
defects and colloidal inclusions, or in frustrated systems. 
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V. CONCLUSION 



Topological defects in directional media are interesting both from the topological and 
geometrical point of view. The physical system cannot undergo a transition between states 
with different topology, which stabilizes the defect states. On the other hand, contributions 
to the free energy that result from spatial variations of the director depend solely on the 
geometry of the field. Constraints conditioned by the free energy are physically just as 
important as topological, which calls for a deeper analysis. 

Instead of treating the nematic director field as a separate entity from the coordinate 
space, we observed it relatively to a chosen surface. Our geometric interpretation of topo- 
logical charge enables directly relating the geometric appearance of the director field to its 
topology. With the texture decomposition, the topological charge can be found simply by 
finding the loops where the director is tangential to the surface and counting them, which 
is a useful representation also for the educational purposes. This manner of calculating the 
topological charge works also for chiral nematics and other media with unit vector fields, 



such as Heisenberg ferro magnets llj and antiferromagnets 23]. 

We show that in nematics, conditions that energetically stabilize the structures are closely 
related to the conservation of surface textures. As the topological charge is ambiguously 
defined in this case, we exploit the geometric nature of the texture decomposition and 
introduce the defect rank, which measures the shape and structure of the director field 
instead of only the topological properties. The defect rank closely matches the intuitive 
understanding of defects and mathematically quantifies the visual clues that are usually 
used to identify the defects by eye. The even or odd parity of the defect rank also matches 
the parity of the topological charge, whenever the defect rank is well-defined. Especially in 
interpretation of the experimental data and description of nematic textures, the defect rank 
should be used instead of the topological charge. The latter should not even be assigned a 
sign when used in a description of a single defect structure, while the defect rank spans an 
entire integer range without ambiguities. This contrasts with the ferromagnetic case, which 
is fully described by the topological charge and does not require description with the defect 
rank. 

Furthermore, the description of defects using surface textures can be applied to systems 



with complex boundary conditions, like nematic colloids 

14 



24 



25|, faceted particles [26], 



nematic droplets [27] and nematics in porous media For a nematic in a homeotropic 
micro- channel network, counting of patches corresponds to the formalism of counting escape 
directions at channel crossings, devised by Serra et al. [3]. The confinement borders of 
the nematic can be viewed as an enclosing surface and the decomposition directly related 
to the boundary conditions enforced by the surface treatment. Colloidal particles with 
heterogeneous surface anchoring also imply a certain texture on its surface 28j. Knowing the 
meaning of the texture features can simplify the design of particles with desired topological 
properties and aid in interpretation of the results. 
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